Incorporating fossils into the joint inference of phylogeny and biogeography of the tree fern order Cyatheales.

Present-day geographic and phylogenetic patterns often reflect the geological and climatic history of the planet. Neontological distribution data are often sufficient to unravel a lineage's biogeographic history, yet ancestral range inferences can be at odds with fossil evidence. Here, I use the fossilized birth-death process and the dispersal-extinction cladogenesis model to jointly infer the dated phylogeny and range evolution of the tree fern order Cyatheales. I use data for 101 fossil and 442 extant tree ferns to reconstruct the biogeographic history of the group over the last 220 million years. Fossil-aware reconstructions evince a prolonged occupancy of Laurasia over the Triassic-Cretaceous by Cyathealean tree ferns, which is evident in the fossil record but hidden from analyses relying on neontological data alone. Nonetheless, fossil-aware reconstructions are affected by uncertainty in fossils' phylogenetic placement, taphonomic biases, and specimen sampling and are sensitive to interpretation of paleodistributions and how these are scored. The present results highlight the need and challenges of incorporating fossils into joint inferences of phylogeny and biogeography to improve the reliability of ancestral geographic range estimation.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns.

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

'Out of Africa' origin of the pantropical staghorn fern genus Platycerium (Polypodiaceae) supported by plastid phylogenomics and biogeographical analysis.

BACKGROUND AND AIMS: The staghorn fern genus Platycerium is one of the most commonly grown ornamental ferns, and it evolved to occupy a typical pantropical intercontinental disjunction. However, species-level relationships in the genus have not been well resolved, and the spatiotemporal evolutionary history of the genus also needs to be explored. METHODS: Plastomes of all the 18 Platycerium species were newly sequenced. Using plastome data, we reconstructed the phylogenetic relationships among Polypodiaceae members with a focus on Platycerium species, and further conducted molecular dating and biogeographical analyses of the genus. KEY RESULTS: The present analyses yielded a robustly supported phylogenetic hypothesis of Platycerium. Molecular dating results showed that Platycerium split from its sister genus Hovenkampia ~35.2 million years ago (Ma) near the Eocene-Oligocene boundary and began to diverge ~26.3 Ma during the late Oligocene, while multiple speciation events within Platycerium occurred during the middle to late Miocene. Biogeographical analysis suggested that Platycerium originated in tropical Africa and then dispersed eastward to southeast Asia-Australasia and westward to neotropical areas. CONCLUSIONS: Our analyses using a plastid phylogenomic approach improved our understanding of the species-level relationships within Platycerium. The global climate changes of both the Late Oligocene Warming and the cooling following the mid-Miocene Climate Optimum may have promoted the speciation of Platycerium, and transoceanic long-distance dispersal is the most plausible explanation for the pantropical distribution of the genus today. Our study investigating the biogeographical history of Platycerium provides a case study not only for the formation of the pantropical intercontinental disjunction of this fern genus but also the 'out of Africa' origin of plant lineages.

Differential expression and co-localization of transcriptional factors during callus transition to differentiation for shoot organogenesis in the water fern Ceratopteris richardii.

BACKGROUND AND AIMS: In flowering plants, regeneration can be achieved by a variety of approaches, and different sets of transcriptional factors are involved in these processes. However, regeneration in taxa other than flowering plants remains a mystery. Ceratopteris richardii is a representative fern capable of both direct and indirect organogenesis, and we aimed to investigate the genetics underlying the transition from callus proliferation to differentiation. METHODS: Morphological and histological analyses were used to determine the type of regeneration involved. RNA sequencing and differential gene expression were used to investigate how the callus switches from proliferation to differentiation. Phylogenetic analysis and RNA in situ hybridization were used to understand whether transcriptional factors are involved in this transition. KEY RESULTS: The callus formed on nascent leaves and subsequently developed the shoot pro-meristem and shoot meristem, thus completing indirect de novo shoot organogenesis in C. richardii. Genes were differentially expressed during the callus transition from proliferation to differentiation, indicating a role for photosynthesis, stimulus response and transmembrane signalling in this transition and the involvement of almost all cell layers that make up the callus. Transcriptional factors were either downregulated or upregulated, which were generally in many-to-many orthology with genes known to be involved in callus development in flowering plants, suggesting that the genetics of fern callus development are both conserved and divergent. Among them, an STM-like, a PLT-like and an ethylene- and salt-inducible ERF gene3-like gene were expressed simultaneously in the vasculature but not in the other parts of the callus, indicating that the vasculature played a role in the callus transition from proliferation to differentiation. CONCLUSIONS: Indirect de novo shoot organogenesis could occur in ferns, and the callus transition from proliferation to differentiation required physiological changes, differential expression of transcriptional factors and involvement of the vasculature.

Unreduced spore formation in a spontaneous chimeric pinnule in an artificially produced haploid Anisocampium niponicum (Athyriaceae, Polypodiales).

Haploid sporophytes of Anisocampium niponicum with 2n = 40, were produced artificially by induced apogamy in vitro. They were subsequently transplanted into pots and two of them have been cultivated for the investigation of sporogenesis and/or production of chimera for more than 20 years. Haploid A. niponicum is sterile, but an abnormal chimeric pinnule that developed spontaneously in a single frond produced sporangia with spores. Each sporangium bore approximately 32 spores that were almost uniform in size. Sowing of these spores resulted in 50 gametophytes. Of 20 gametophytes cultured individually, five produced sporophytes apogamously after eight months. Both the gametophytes and subsequent apogamous sporophytes showed a chromosome number of 2n = 40. Our study demonstrates that a haploid sporophyte offspring can be produced from a haploid mother sporophyte via haploid spores. Since asexual reproduction is a prominent evolutionary process in ferns, the reproduction of a haploid A. niponicum sporophyte by unreduced spore formation might help to elucidate how apogamous ferns occur and evolve.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns.

BACKGROUND: Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events. RESULTS: In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species. CONCLUSIONS: Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Fern proteins show promise against crop pests.

These and other ancient plants may provide alternatives to chemical insecticides.

The evolutionary innovation of root suberin lamellae contributed to the rise of seed plants.

Seed plants overtook ferns to become the dominant plant group during the late Carboniferous, a period in which the climate became colder and dryer1,2. However, the specific innovations driving the success of seed plants are not clear. Here we report that the appearance of suberin lamellae (SL) contributed to the rise of seed plants. We show that the Casparian strip and SL vascular barriers evolved at different times, with the former originating in the most recent common ancestor (MRCA) of vascular plants and the latter in the MRCA of seed plants. Our results further suggest that most of the genes required for suberin formation arose through gene duplication in the MRCA of seed plants. We show that the appearance of the SL in the MRCA of seed plants enhanced drought tolerance through preventing water loss from the stele. We hypothesize that SL provide a decisive selective advantage over ferns in arid environments, resulting in the decline of ferns and the rise of gymnosperms. This study provides insights into the evolutionary success of seed plants and has implications for engineering drought-tolerant crops or fern varieties.

Leaf structure and water relations of an allotetraploid Mediterranean fern and its diploid parents.

Allopolyploidy is a common speciation mechanism in plants; however, its physiological and ecological consequences in niche partitioning have been scarcely studied. In this sense, leaf traits are good proxies to study the adaptive capacity of allopolyploids and diploid parents to their respective environmental conditions. In the present work, leaf water relations (assessed through pressure-volume curves) and structural and anatomical traits of the allotetraploid fern Oeosporangium tinaei and its diploid parents, Oeosporangium hispanicum and Oeosporangium pteridioides, were studied under controlled conditions in response to a water stress (WS) cycle. O. hispanicum showed the lowest osmotic potential at turgor loss point (πtlp ) and leaf capacitance, together with higher leaf mass per area (LMA), leaf thickness (LT), leaf density (LD), and leaf dry matter content (LDMC), whereas O. pteridioides presented the opposite set of traits (high πtlp and capacitance, and low LMA, LT, LD, and LDMC). O. tinaei showed an intermediate position for most of the studied traits. The responsiveness (osmotic and elastic adjustments) to WS was low, although most of the traits explained the segregation of the three species across a range of drought tolerance according to the rank: O. hispanicum > O. tinaei > O. pteridioides. These trait differences may underlie the niche segregation among coexisting populations of the three species in the Mediterranean basin.

Expression and Functional Studies of Leaf, Floral, and Fruit Developmental Genes in Non-model Species.

Researchers working on evolutionary developmental plant biology are inclined to choose non-model taxa to address how specific features have been acquired during ontogeny and fixed during phylogeny. In this chapter we describe methods to extract RNA, to assemble de-novo transcriptomes, to isolate orthologous genes within gene families, and to evaluate expression and function of target genes. We have successfully optimized these protocols for non-model plant species including ferns, gymnosperms, and a large assortment of angiosperms. In the latter, we have ranged a large number of families including Aristolochiaceae, Apodanthaceae, Chloranthaceae, Orchidaceae, Papaveraceae, Rubiaceae, Solanaceae, and Tropaeolaceae.

A Ceratopteris EXCESS MICROSPOROCYTES1 suppresses reproductive transition in the fern vegetative leaves.

Land plant sexual reproduction involves the transition of cells from somatic to reproductive identity during post-embryonic development. In Arabidopsis, the leucine-rich repeat receptor-like kinase EXCESS MICROSPOROCYTES1 (EXS/EMS1) restricts the number of sporogenous cells during the transition from diploid tissue to haploid spore production by promoting the formation of the tapetum cell layer within the anther. Although all land plants studied contain EMS1 genes, its function is unknown beyond a few angiosperms. In the model fern Ceratopteris (Ceratopteris richardii), we discovered an EMS1 homolog (CrEMS1) that functions to suppress formation of reproductive structures on vegetative leaves of the fern sporophyte, a role not found in angiosperms. Suppression of CrEMS1 by RNAi did not affect sporogenesis on reproductive leaves but did affect antheridium production of the fern gametophyte. Expression patterns of CrEMS1 across developmental stages suggest threshold levels of CrEMS1 control the specification of reproductive organs during both generations of the fern. Additional EMS1 homologs present in the fern genome suggest a dynamic role of EMS1 receptors in the evolution of reproductive development in vascular plants.

Proteome and Interactome Linked to Metabolism, Genetic Information Processing, and Abiotic Stress in Gametophytes of Two Woodferns.

Ferns and lycophytes have received scant molecular attention in comparison to angiosperms. The advent of high-throughput technologies allowed an advance towards a greater knowledge of their elusive genomes. In this work, proteomic analyses of heart-shaped gametophytes of two ferns were performed: the apomictic Dryopteris affinis ssp. affinis and its sexual relative Dryopteris oreades. In total, a set of 218 proteins shared by these two gametophytes were analyzed using the STRING database, and their proteome associated with metabolism, genetic information processing, and responses to abiotic stress is discussed. Specifically, we report proteins involved in the metabolism of carbohydrates, lipids, and nucleotides, the biosynthesis of amino acids and secondary compounds, energy, oxide-reduction, transcription, translation, protein folding, sorting and degradation, and responses to abiotic stresses. The interactome of this set of proteins represents a total network composed of 218 nodes and 1792 interactions, obtained mostly from databases and text mining. The interactions among the identified proteins of the ferns D. affinis and D. oreades, together with the description of their biological functions, might contribute to a better understanding of the function and development of ferns as well as fill knowledge gaps in plant evolution.

Non-seed plants are emerging gene sources for agriculture and insect control proteins.

The non-seed plants (e.g., charophyte algae, bryophytes, and ferns) have multiple human uses, but their contributions to agriculture and research have lagged behind seed plants. While sharing broadly conserved biology with seed plants and the major crops, non-seed plants sometimes possess alternative molecular and physiological adaptations. These adaptations may guide crop improvements. One such area is the presence of multiple classes of insecticidal proteins found in non-seed plant genomes which are either absent or widely diverged in seed plants. There are documented uses of non-seed plants, and ferns for example have been used in human diets. Among the occasional identifiable toxins or antinutritive components present in non-seed plants, none include these insecticidal proteins. Apart from these discrete risk factors which can be addressed in the safety assessment, there should be no general safety concern about sourcing genes from non-seed plant species.

The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation.

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Splitting one species into 22: an unusual tripling of molecular, morphological, and geographical differentiation in the fern family Didymochlaenaceae (Polypodiales).

The pantropical fern genus Didymochlaena (Didymochlaenaceae) has long been considered to contain one species only. Recent studies have resolved this genus/family as either sister to the rest of eupolypods I or as the second branching lineage of eupolypods I, and have shown that this genus is not monospecific, but the exact species diversity is unknown. In this study, a new phylogeny is reconstructed based on an expanded taxon sampling and six molecular markers. Our major results include: (i) Didymochlaena is moderately or weakly supported as sister to the rest of eupolypods I, highlighting the difficulty in resolving the relationships of this important fern lineage in the polypods; (ii) species in Didymochlaena are resolved into a New World clade and an Old World clade, and the latter further into an African clade and an Asian-Pacific clade; (iii) an unusual tripling of molecular, morphological and geographical differentiation in Didymochlaena is detected, suggesting single vicariance or dispersal events in individual regions and no evidence for reversals at all, followed by allopatric speciation at more or less homogeneous rates; (iv) evolution of 18 morphological characters is inferred and two morphological synapomorphies defining the family are recognized-the elliptical sori and fewer than 10 sori per pinnule, the latter never having been suggested before; (v) based on morphological and molecular variation, 22 species in the genus are recognized contrasting with earlier estimates of between one and a few; and (vi) our biogeographical analysis suggests an origin for Didymochlaena in the latest Jurassic-earliest Cretaceous and the initial diversification of the extant lineages in the Miocene-all but one species diverged from their sisters within the last 27 Myr, in most cases associated with allopatric speciation owing to geologic and climatic events, or dispersal.

Systematics and biogeography of the Old World fern genus Antrophyum.

Antrophyum is one of the largest genera of vittarioid ferns (Pteridaceae) and is most diverse in tropical Asia and the Pacific Islands, but also occurs in temperate Asia, Australia, tropical Africa and the Malagasy region. The only monographic study of Antrophyum was published more than a century ago and a modern assessment of its diversity is lacking. Here, we reconstructed a comprehensively sampled and robustly supported phylogeny for the genus based on four chloroplast markers using Bayesian inference, maximum likelihood and maximum parsimony analyses. We then explored the evolution of the genus from the perspectives of morphology, systematics and historical biogeography. We investigated nine critical morphological characters using a morphometric approach and reconstructed their evolution on the phylogeny. We describe four new species and provide new insight into species delimitation. We currently recognize 34 species for the genus and provide a key to identify them. The results of biogeographical analysis suggest that the distribution of extant species is largely shaped by both ancient and recent dispersal events.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns.

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

Ecological distribution and genetic diversity of Azolla in Uganda.

BACKGROUND: Azolla is an important aquatic fern whose agronomic potential has not been fully exploited in Uganda. This study aimed at determining the genetic variation in the Azolla species existing in Uganda and the factors influencing their distribution in the different agro-ecological zones of Uganda. Molecular characterization was preferred in this study because of its efficiency in detecting variations among closely related species. RESULTS: Four species of Azolla were identified in Uganda with 100, 93.36, 99.22 and 99.39% sequence identities to the reference database sequences of; Azolla mexicana, Azolla microphylla, Azolla filiculoides and Azolla cristata, respectively. These different species were distributed in four out of the ten agro-ecological zones of Uganda which are situated in close vicinity to large water masses. The principal component analysis (PCA) results revealed that maximum rainfall and altitude significantly accounted for the variations in the distribution of Azolla with factor loadings of 0.921 and 0.922, respectively. CONCLUSION: Massive destruction coupled with prolonged disturbance of Azolla's habitat negatively affected its growth, survival and distribution in the country. Therefore, there is a need to develop standard methods that can preserve the various species of Azolla, so as to salvage them for future use, research and reference.

Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis.

Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.

Fern cell walls and the evolution of arabinogalactan proteins in streptophytes.

Significant changes have occurred in plant cell wall composition during evolution and diversification of tracheophytes. As the sister lineage to seed plants, knowledge on the cell wall of ferns is key to track evolutionary changes across tracheophytes and to understand seed plant-specific evolutionary innovations. Fern cell wall composition is not fully understood, including limited knowledge of glycoproteins such as the fern arabinogalactan proteins (AGPs). Here, we characterize the AGPs from the leptosporangiate fern genera Azolla, Salvinia, and Ceratopteris. The carbohydrate moiety of seed plant AGPs consists of a galactan backbone including mainly 1,3- and 1,3,6-linked pyranosidic galactose, which is conserved across the investigated fern AGPs. Yet, unlike AGPs of angiosperms, those of ferns contained the unusual sugar 3-O-methylrhamnose. Besides terminal furanosidic arabinose, Ara (Araf), the main linkage type of Araf in the ferns was 1,2-linked Araf, whereas in seed plants 1,5-linked Araf is often dominating. Antibodies directed against carbohydrate epitopes of AGPs supported the structural differences between AGPs of ferns and seed plants. Comparison of AGP linkage types across the streptophyte lineage showed that angiosperms have rather conserved monosaccharide linkage types; by contrast bryophytes, ferns, and gymnosperms showed more variability. Phylogenetic analyses of glycosyltransferases involved in AGP biosynthesis and bioinformatic search for AGP protein backbones revealed a versatile genetic toolkit for AGP complexity in ferns. Our data reveal important differences across AGP diversity of which the functional significance is unknown. This diversity sheds light on the evolution of the hallmark feature of tracheophytes: their elaborate cell walls.